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The exceedancc point process approach of 
Hsing Cl al. is extended to multivariate 
stationary sequences and some weak con- 
vergence results are obtained. It is well 
known that under general mixing assump- 
tions, high level exceedanccs typically 
have a limiting Compound Pbisson structure 
where multiple events arc caused by the 
clustering of cxccedances. In this paper we 
explore (a) the precise effect of such 
clustering on the limit, and (b) the relation- 
ship between point process convergence 
and the limiting behavior of maxima. Fol- 



lowing this, the notion of multivariate ex- 
tremal index vs. introduced which is shown 
to have properties analogous to its uni- 
variate counterpart. Two examples of 
bivariate moving average sequences are 
presented for which the extremal index is 
calculated in some special cases. 
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1. Introduction 



Extreme value theory for multivariate iid sequences 
has been studied for quite some time now but attention 
to the dependent case has been relatively recent. For 
univariate sequences it is known that local dependence 
causes extreme values to occur in clusters, which in turn 
results in a stochastically smaller distribution for the 
maximum than if the observations were independent. We 
begin with a brief review of these results, which we 
shall later extend to the multivariate case. 

Let {£, } be a univariate stationary sequence. Write 
Af„=max {£,,...,£„} and for t>0, let {u„(t)} denote :■ 
sequence satisfying lin%,^/i/ J {£i>h 1 ,(t)}=t. Under 
quite general mixing assumptions there exist constants 
0=£f)'<t?"£l such that 

lini sup ^{A/^h^tJHc "' and 

liminfP<M^w„(T)}-e-*" T 



for all t. (See Ref. [1], although the idea actually dates 
back to Refs. [2-4].) Thus if P {M „<u„(t )} converges 
for some t , then d'=ff"(=6, say) and hence 
lim„^ f , {A/*<K„(T)}=e"*' for all r>0. The common 
value f? is then called the extremal index of {£, }. We 
shall assume 9 to be positive whenever it exists, since the 
case 6=0 corresponds to a degenerate limiting distribu- 
tion for M„ . Note that 6=\ for iid sequences. Let {£, } be 
an iid sequence with £i— ** £i , called the associated iid 
sequence, and write A/„=max{^i,...,^„}. If {£,} has ex- 
tremal index 6 and lim„^ P{M n ^v„(t)}=H(t) for a 
suitable family of normalizing constants {v„(r)}, then it 
follows (upon identifying e"* r with //(/)) that lim,^ 
P{M H <v a (t)}=H(t) where 

H{t)=H{t)\ (I) 

The extremal index is thus a measure of the effect of 
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dependence on the limiting distribution of A/„. The 
stochastically smaller limiting distribution of A/„ is in 
fact a direct result of the clustering of extremes, as 
explained below. See Ref. [5] for details. 

For fixed 7>0 let the e.xceedance point process 
N n ^Nl n be defined by 

JV„(£)-£ /(^mjWw. BC[0,\], 
i-i 

where l A denotes the indicator function of the event A . 
Then for a broad class of weakly dependent sequences, 
the limit in distribution of rV,, if it exists, is a Compound 
Poisson process with intensity 0r and multiplicity distri- 
bution 7ron {1,2,...}. 'llie Poisson events may in fact be 
regarded as the positions of "cxccedancc clusters" 
while the multiplicities correspond to cluster sizes. 
More explicitly, one may divide the n obsei"vations into 
Jt„ blocks of roughly equal size and regard exceedances 
within each block as forming a single "cluster", so that 
the cluster sizes arc given by N m (J,). i=l,...Jc„, where 
W^M'f , j£J- Hot a suitable choice of Jt„ depending 
on the mixing rate of {£, >, one then has 

l\m P {N„(Ji>j\ NJJ,)>0}*nU). J21, 

and 

lim P'"{/V„(./ l )-OH™ P{N.l*),n=0> 

-lim P {*/„=£ w„(r)}=c -*, 

so that in particular, lmi„_„, k„ P {N a Ui)>O}~0r, Hence 

liin£A r „l().IHimJt./i"A',v7 l ; 

-flrlim B(JV,y,)|/V.y,)>f/>. 

Ji->=e 

while on the other hand, lini„_^ A'A r ,l(),l|=lini„_ r 
nP{(\>u„(T)}~T. The cluster size distribution and the 
extremal index are therefore related by 



li m E(AW,)|A r ,(/,)>OH/fl. 



(2) 



Now let {in=(£nu—.Cd), n£E2Z} be a multivariate sta- 
tionary sequence where d>\ is a fixed integer, and 
write M„-(M„ U -Mn,i) where A/„/=iiiax{6,-, ...,£,,}, 
j—]„..,d. The study of multivariate extremes began in 
the early 1950s, focusing mainly on the limiting behav- 
ior of hi, under a linear normalization, when the obser- 
vations arc iid. The resulting class of limiting distribu- 
tions was characterized in Ref. [6] and domains of 
attraction criteria were given in Ref. [7J. See also Ref 



[8J. Chapter 5, for an account of the literature surround- 
ing this theory. For stationary sequences satisfying a 
general mixing assumption, it is known (see Rcfs. 19, 
10], and Theorem 1.1 below) that the class of limiting 
distributions of M„ is the same as for iid sequences. In 
this paper we explore the precise effect of dependence 
on the limiting distribution by extending the univariate 
theory described above to the multivariate case. Essen- 
tially, this involves studying the interrelationship be- 
tween the two dependence structures present, one due to 
dependence over time and the other due to the depen- 
dence between the various components of the multivari- 
ate observations. The ideas become most transparent 
when presented in terms of so-called dependence func- 
tions IHj. Here we adopt the slightly modified definition 
found in Ref. [9]. A distribution function D on [0,1V is 
called a dependence function if Dj(Dj{u))=D t (u), 
w£l<),l |, j=l,...,d. where the subscript y signifies they'th 
marginal. The dependence function of a distribution F 
on IR J is defined by D f (K)=P{F,(X,)< 

hi F^Xj^u,,}, «=(«i «rf)£ 10J] J , where 

(X u ...J(j) is a random vector with distribution F. More 
generally, any dependence function satisfying 
F(x)*D(F\(x,),.,.J'',tiXj)) could be defined to be a de- 
pendence function of F, although the present choice is 
a natural one. 

Write r={0,l)A{I> where l=(l l)eiR J , and for 

t=(r 1 ,....^)er, let f,(l(Hr,i(/i),.,v^)) where i'„ ; (/j) 
satisfies lim„_,« «/'{£i > >i' lu (;,)}=— log( ; . Let W„ denote 
the distribution function of M n (i.e., H„(x)*-P {M,^x }), 
with marginals //„,-, >*!,...,(/". Then (see Refs. fS, II]). 



H,li>,(t)H*Hlt) 



0) 



if and only if 



The limiting behavior of M„ can therefore be separated 
into two parts, one pertaining to the convergence of the 
marginals (a univariate problem) and the other to the 
convergence of the dependence functions. Here we fo- 
cus attention exclusively on the latter. It should be noted 
that the choice of normalising constants does not affect 
the dependence function of the limit distribution H , hut 
only alters the marginals (sec Ref. (91, Lemma 3.2). 
Since our main interest is in the dependence function, 
the present choice of normalising constants is appropri- 
ate in view of the fact that it results in Uniform[0, 1]' 
marginals for the limit distribution when {£, > is iid, so 
that in particular Dn=fl . According to Theorem 3.3 of 
Ref. [91, the class of all possible limits H in Eq. (3) (for 
iid {£,}) Ls precisely the class of extreme dependence 
functions, that is those that satisfy 
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D"(t)-0(r,\...,rj) 



(4) 



A^)-2Wa>S„ aeio,i). 



(5) 



for each n>l andt=(f j d )£[0, 1 Y. Theorem 1.1 below 

shows that the same is true also if {£, } is a stationary 
sequence satisfying the following mixing condition. 
For ter. let 

HH*kWh<r<(&>v<jlif)) ■ ksisl,/-l d} 

and for 1 ^/^«— I , define 

a„j-Sup{\p(Ar)B)-P(A)P(B)\ : A€9§f(»„(t)), 

BeSC(v„(t)), 1 <jt<*-t-/^rc }. 

The mixing condition A(*'„(t}) is then said to hold if 
a, Jn — »0 for some sequence {/„ } satisfying /./«— >0. This 
is the multivariate version of the mixing condition used 
in Ref. [51 and is slightly stronger than the D(u„) condi- 
tion in Ref. [91- Henceforth {£„ } will be assumed to 
satisfy A(v„(t)), for some or all t, as required. 

Theorem 1.1. Let {&} satisfy A(v ■„(!)) for all tEf and 
suppose that P{M n ^v„{\))— >"' H(t), non-degenerate. 
Then D„ is an extreme dependence function and hence, 
in particular, #(t<)=ff'{t) for each tE[Ojf and e>0 
(where f-(tf,...,t^). 

Proof: The first part is an immediate consequence of 
Theorem 4.2 of Ref. [9] while the second part follows 
from the definition of extreme dependence functions 
upon noting that (by the univariate theory described 
above), the marginals of H are of the form Hj{tj)=tf' 
where fy is the extremal index of {£„j}, the _/th -compo- 
nent sequence of {£, }. 

In the next section we apply the cxcccdance point 
process approach to multivariate extremes and obtain 
some weak convergence results. The multivariate ex- 
tremal index is then defined (in Sec. 3). based on the 
multivariate analogue of Eq. (1). It is seen to be a func- 
tion of only d— 1 variables and its properties naturally 
extend those of the univariate extremal index. Finally in 
Sec. 4 we consider two examples of bivariatc moving 
average sequences for which the computation of the 
extremal index is demonstrated. 

2. Exceedance Point Processes 

Fix rG7"and let ^ifh^-.p^ '=' ,...jt.j-l,...i, and put 
S,=(5,i,...,8i/). The multivariate exceedance point pro- 
cess N^N™ is then defined by 



Assume that {f„} satisfies A(v fl (t)). If also M,-/ Mi. 
then it may be shown (as in the univariate case) that the 
limit N n is a point process on [0,1] which is of Com- 
pound Poisson type. More precisely, the Laplace Trans- 
form of M) is given by 



-log E cxp{- 2 \ fj dN 0i )-v\ J 
/-i J |o.i| •'jeio.jj-'j 



e*f 



(l-exp{-2y/;C*)»d7r(y)dv. 



(6) 



Here N j denotes the jth-componcnt of No, v is a 
positive constant, it is a probability distribution on 
2Lf =-{0,1,2,... y\{0} andj^'s are non-negative functions 
on [0,11. 

Let {k„ } be any sequence of positive integers satisfy- 
ing 

* B ->^, kJJn-^0, and k*& hJa ->0 r . ,*$.*-vx. (7) 

Set r„=[n!k„] (the largest integer not exceeding nlk„) 
and put J„i=[0,f„/«[. Define the probability distribution 
77, on Z? by 

*.(y W» {N„ (A,)-y I M, (•/,!)*«}, yesf . 

The following theorem which gives a useful characteri- 
zation of the convergence of N„ is an immediate conse- 
quence of the results in Sec. 5 of Ref. [12]. 

ThkOiu-m 2.1, M„-/ jV if and only if tt„— >"' tr and 
P {Af„^v„(t)}— > e~*. and in that case the Laplace Trans- 
form ofN a is given by Eq. (6). 

Next we consider the iid case in some detail and 
obtain an interesting connection with Theorem 5.3.1 of 
Ref. [8]. 

Proposition 2.2. Let {£,} be iid and for fixed tET let 
N n be defined by Eq. (5). lfN„—>'' N then the multiplic- 
ity distribution it in Eq. (6) is supported on the set 
S-{0.1 }\{0}. 

Proof: Observe that A(v„(t)) is trivially satisfied since 
a„.i=0 so that we may take /„sl and A„«n. Then 
w.(jW* {<5,=y 1 fi^v„(t)}, fE2Zt. which is clearly sup- 
ported on S. The result is now immediate since 5 is a 
closed set and ir„— ** tt by Theorem 2.1. 

Making the dependence on t explicit, wc now write 
/v><", MWVj", v~v w and ttW 1 . In addition we shall 
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require the following notation from Ref. [8], For 
l<it<rf, let j(A)=0 J*) denote a vector with integer- 
valued components lsc/i<&<— <jt£d, and for 
x=(Xi,...jc rf )£lR'' write xw-{x h ,..*Xj t ). Define the "sur- 
vival function" 

G (x)-P {£, ,>*„ .. .,ftj>Jr rf } 

and write G K i)(x)=P<^,>J ;i ,...,f[ A >Jf ;i }. For each j(k), 
let yj (4) denote the element in 5= {0,1 }"\{0} whose jth 
component equals 1 if and only if/*/,- for some i- 1 ,...,k. 
(This defines a natural 1-1 correspondence between 5 
and the j(Jt)'s.) 

Theorem 2.3. Let {£, } be iid. Then NS^-i'N^for some 
fixed tGT" if and only if 



lim RGjmfoan-AjtiXtVe" 



(8) 



/or eoc/t j(Jt). l<A-<i/. In rAo/ awe A^" /iot Laplace 
Transform given by Eq. (6) with 



tfrtd wif/r 7T ( ° determined by the relations 



(9) 



WtW 2 **"&)■ CO) 

Proof: Write S^O-St^,^*/?^!), so that P{S,£ 
*»(t)>«XJ.,(-l) l * , Si<v,<t)). If Eq. (8) holds for each *, 
then 

4 

\imnP<&£v n (t)}=2(-l) M 2 Awjtt)-^', (11) 

and hence lim 1M «/ > {Af lt Sv 11 (t)}-e~ l/ ' > . 
Next observe that for each j(A'). isJc^d, 

G x dv-.m>- 2 nfajy-p <&**$» E ^(y),<i2) 

where 7r„(y)-/'{6i-y I f,*v„(t)}. Moreover, this rela- 
tionship is invertiblc in the sense that each of the prob- 
abilities -n-„(y), y&S, can be expressed as a linear com- 
bination of the Gj(n(v„(t))*s. Therefore by Eq. (8). 
lim„_ rt 7r J ,(y)=-7r ( "(y) (say) exists and satisfies Eq. (10). 
Hence by Theorem 2.1 A?'-*'' Mf where Aft' has the 
specified parameters. Conversely if A'"'— » J A^" then tr n 
and P{A/„<v„(t)} converge (by Theorem 2.1 again), 
and hence Eq. (8) follows by virtue of Eqs. (II) and 
(12). D 



Corollary 2.4. Let {&} be iid. Then N™-*' Aft 1 for 
each ter if and only if P {M„=£v n (t)}-» H " H(t). More- 
over H and <v l,1 ,7r" ! } ie j- determine each other. 

Prook: (Sketch) The first part follows from Theorem 
2.3 above and Theorem 5.3.1 of Ref. [8] which states 
that PW n ^v„(t)}^ v H(t) if and only if Eq. (8) holds 
for each tG7". Note that //(i)^"*"' so that H and the 
v m 's can be obtained from each other. Also the ir m 's 
can be obtained from the ^"'s by first inverting Eq. (9) 
to get the ftj(i){t)'s and then inverting Eq. (10). (The 
inversion of Eq. (9) is carried out inductively using the 
fact that the weak convergence of rV„(v„(t)) implies that 
of all lower dimensional marginals.) □ 

Analogous results for the dependent case take on a 
slightly different form. Let {£, } be a stationary se- 
quence satisfying A(v J ,(t)) for each tG7\ As before let 
r»=[rt/A,J where {k„} is any sequence satisfying Eq. (7). 
and define 

Gr^(v„ (t))=P {M,j ,>V. j,(t ,),.~Mrj t > I'™ (f» ) >- 

Theorem 2.5. Lei {£„ }be a stationary sequence satis- 
fying &(v„{t))for each ter Then P {M„^,,(t)}-r H (t) 
if and only if 

limk„G, M ,(v n (t))~h m (t)<™ 

n-*oc 

for each j(Jt), 1 ^k^d and tGT, and in that case 

4 

//(t)=exp{2<-l)' X **#»■ 

Proof: Observe that the mixing condition A(v„(t)) im- 
plies that {f„j, () } satisfies A(v„j< () (t)) for each j(Jfc) (with 
obvious notation). Hence it may be shown as in the 
univariate case (see Lemma 2.1 of Ref. (11) that 

P {M„ j(() <v^ i ,(l)}-P t * <M ri ji»*».j*(t))-»0, (13) 

for each j(A). The result may therefore be proved in 
exactly the same way as Theorem 5.3.1 of Ref. [8]. □ 

Remark: Under the hypothesis of Theorem 2.5, if 
Atf'-y* Aft', tGT", with parameters i> m and tt" 1 then 
P{M„<i',(t)}-^" //(tj-c--* 1 , as in the iid case. However 
it is not possible in general to recover the ir' n 's from H 
since the clustering of exceedances may cause the sup- 
port of tt ( " to extend beyond S. References [9, 10] give 
sufficient conditions (analogous to the D'(u„) condition 
of Ref. [13]) under which clustering docs not occur, so 
that Corollary 2.4 can be extended to stationary se- 
quences satisfying this condition. 
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A distribution function F on JR rf is said 10 be indepen- 
dent if F(x)-ir H Fj(.Xj), x&R d . Jf {£,} is iid awl 
P {M„^v„(t)}-> v /7(f). then ii follows from CoioHary 
5.3.1 of Ref. [8] thai H is independent if and only if die 
marginals of H are pairwise independent. The analogous 
result for the dependent case is stated below. The pi oof 
(which is omitted) is essentially the sanie as for the iui 
case, but uses Theorem 2.5 instead of Theorem 5.3. t of 
Rcf. [8]. 



3. The Multivariate Extremal Index 

i ei {£„ I he a Baiiomiry sequence and {£„ } the associ- 
ated nd secjueive. Suppose that P{M.^v R (t)}-* w H(t) 
and f'{M^v t {t)} -»"' H(t). The multivariate extremal 
indsx ot {£,} is (hen defined by the relation 
H{\)-ll M (i) (see Eq. (1)), or more explicitly 

auHfj* //(t>/iog //(t), ter. 



CorOliary 2.6. Z-e/ {£, } fee a stationary sequence ult- 
isfying A(v„(t)) for each tE!T and suppose ih.it 
P {W J ,<v,(t)}->" //(t). T/hm tf ii independent if and 
only if k„P{M,j>v„j{tj) t M r j>v a j(t,)}-^0 for each 
lSj<l3Ul. VET, i.e., if and only if k„G flMk) (vM) *Ojar 
each j(2) ant/ each IE.T. 

It is shown in Ref. [14] that H is independent it 
HityUj.iHjitj) for some tG(0,)) J . Although the result 
in [14] only stated for iid sequences under a linear nor- 
malization, the proof essentially rests on the defining 
property of extreme dependence functions, namely Eq. 
(1). Consequently the result extends lo the present more 
general situation allowing dependence and non-lineai 
normalizations. Corollary 2.6 can therefore be improved 
as follows. 

Corollary 2.7. Let {§, } be as in Corollary 2.6 and 
suppose that P {M H <v„(t)}-* w H(t). Then the fallowing 
are equivalent: 
(i) H is independent, 
(ii) H(V r ~n% l H,(t i )for some te(0,l)*, 
(Hi) *„G roKtl (v,(t))^0/or each }(2).for some t£(0,l) J . 

It should be noted that Refs. [9, 10] give somo inter- 
esting sufficient conditions for H to be independent 
when {£,} is a stationary sequence. A natural question 
to ask in the present context is whether H is independent 
whenever H is. Proposition 3.4 gives a necessary and 
sufficient condition for this in terms of the extremal 
index, but the answer in general is negative and a coun ■ 
ter-examplc can be found in [10]. It seems more plausi- 
ble that the converse may be tme, i.e., that H is indepen- 
dent whenever H is. In fact however, this loo is not the 
case, as shown by an interesting counter-example in 
[15]. 

We conclude this section by stating a result which 
extends Theorem 5.1 of [5| and is proved similarly. 

Thhorf.m 2.8. Let {£, } be a stationary sequence satis- 
fying &{v,(t))for each t&T and suppose thatN^Stt"'' 
for some t£T. Then JV<' c, -»' f N^ for each »0 and fur 
thermore, /W/ 1 and tt^W 1 (where ?-{tf ,...$)). 



Obstnc that frit) is well defined since H has Uni- 
fo>m[0,)| marginals and hence, 0<//(t)<l on T. The 
following results describe some basic properties of the 
multivariate extremal index. 

I'ROPOSnxJN 3.1. Assume thai {£,} satisfies A(v„(t)) for 
each MET and has extremal index 9(\). Then 
(i) 0(t)^6(t r ) for each t^T and c>0, and 
(ii) foi each y'-l „..,£/, {£„;} has extremal index 6j=0(t) 

where lEThas all coordinates equal to I except the 

jth. 

(Note that by (i), B{\) is constant along the contours 
l,.-*{l : u'>0>, ttT, and hence Bj in (//) is well-defined.) 

I'sook Recall that (by Theorem 1.1) //(0-// c (t) and 
/ytf')«# ; (t) so that (/) foUows from the definition of the 
extremal index. Next, for t£7" with all coordinates but 
ihe/th equal to 1, P{M K ^v„(t)}-P {M^v^tf)} and 
hence 



1 



im />{««,* i^(iy) Him /»{Jt/.*r.(t)W/(t)-«,(0). 

ITieicfure by Theorem 2.2 of Ref. [1], {£ nJ } has ex- 
tremal index d, (say) so that Hj(tj)=tp. Now 
/l(i)"-ff m Xt) by definition of the extremal index, and for 
the present choice of t this is the same as Hj(tj}=t* l) , 
whence it follows that 6(1)=$ for all such t. □ 

For ter, let JVi" denote the one-dimensional point 
process obtained from iVi" via the map y— *l<j± i from 

{o.i}' to {o.i}, i.c., w«»(fl)-2ti i^j***), Bern. 

Thus A'"' has unit mass at iln if and only if £^v„(t). 
Assume that {£, } satisfies A(»'„(t)) and with ■/„, as in 
Sec. 2. lei 



<"(y)-/> (A?rV.i)-y 1 NH>(J„,)>0}, y^l. 

Pw.*OS11"K'jn 3.2. Assume that {£,} satisfies A(v„(t)) for 
each IE.T and has extremal index $(t). Then 
0(ty^(lim i ,.,JL^ ] y^'iy))-'. 
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Proof: Observe that 



r„/M6£v„<t)} 



k„r a nP{g,^v.(t)> 
n Jt a P{A/ r „Sv,(t)}' 

Now linw, P{M n <v„(t)}-W(t) and lim,^ 
/>{A/„<v»(t)}-/f(t) (by assumption), so that Iim„_^ 
«P{£,3Sv„(t)} — log H(t) and (by Eq. (13)) Vinw 
^/ > {W,,£v„(t)}=-log H(t). Therefore lim,^. l f « t 
y^"(y)=Iog #(t)/log //(t)-l/0(t), as required. D 

Rbmakk: Proposition 3.2 is simply the multivariate ver- 
sion of Eq. (1) and shows how the extremal index is 
related to the clustering of "exceedances." Indeed, ac- 
cording to the present viewpoint, an excccdance occurs 
at time i if £*v„(t), i.e., if £,j>v BJ (/j) for at least one j. 
Thus Propositions 3.! and 3.2 show that while the de- 
gree of clustering may depend on t, it is constant on each 
L,. Note also the connection to Theorem 2.8. 

The next result gives the relation between the depen- 
dence functions of H and H, which is seen to involve the 
extremal index in an intricate manner. 

Proposition 3.3. // {£,} has extremal index 0(t), t£T. 
then 



D H (t* /;<)=D£«>(t), ter, 

where Bj is the extremal index <>/{§,;}, j=\ d. 



(H) 



Proof. By definition of the dependence function, 
D H (tK{//,(X,)</ l ,..AX)<a where (X,,...X) is 
a random vector with distribution H . Therefore, since 



'} * 



D„W tf)=P{X>^t u ...J(^t J }=H(t)-H i "Xt) 



-D,f'(t). 



or required. D 



Rhmarks 

] ,) Note that s=t c (for some oO) if and only if log s^/log 

sj=\og fj/log tj-aj (say), j=\ d—\. Therefore we 

may write L,-L, where a=(fl|,...,aj_i), and hence by 
the remark following Proposition 3.2, 0(t)«*#(a), 
i.e., the extremal index is a function of d— 1 vari- 
ables only. 

2.) By Proposition 3.3, D w (f 1 fl ',...,/ < ?0=^' J (t)-D fl (t*"). 

Also, if t£L, then (/*' ^')eL,. where a*=(«j0|/ 

9j a a -\6 d -il6,i). Thus D H is obtained by translat- 
ing the values of D,;(=ft) on 1, onto L a .. 



3.) While the above results illustrate some of the basic 
properties of the multivariate extremal index, they 
are far from complete. For instance; it would be 
useful to identify the set of all "admissible" &(•) for 
a given ft, that is the set of all #(■) such that D„() 
defined by Eq. (14) is a probability distribution on 
[0,1]''. It would also be of interest to study the prop- 
erties of #(■) when one or both of ft and H are 
independent. In this context we haw the following 
simple result. 

Proposition 3.4. If ft is independent, then If is inde- 
pendent if and only if 

d d 

0(t)-2 6, logfj/2 lo £0 . M some t £(0,1 )", 

In particular, if both ft and H are independent then 8(t) 
is a convex combination of the 9/s. 

Proof: If ft is independent, then H(t)-H(t) m = 
(II^,(j)** J . The conclusion follows immediately from 
Corollary 2.7 (iii) upon taking logarithms and noting 
that if H is independent, then //(t)=n£ 1 //'- 

Thc extremal index can be given the following more 
general formulation. Let /I and /x be the probability 
measures on (0,1)'' corresponding to ft and H , respec- 
tively. Thus for instance, 

where v„(A)-{v«($) : sGA >, AC(OJ)''. We now define 
6(A) via the relationship ix(A)-jl WA, (A), or more di- 
rectly e(.A)=\og/jL(A)/]o«p.(A), for subsets ^C(0,l) J 
such that p.(A)>0 and /i(A)>0. 

Note that 0(t)-0((O,/,)X -XfO,/,,)) for tG7. Thus if 
{0(t) : tET} is known along with either of// or ft, then 
it is possible at least in theory to obtain {0(A) : 
ACZ(0,\) d }. In practice, however, it may not be possible 
to obtain 0(A) in a tractable form, but frequently one is 
only interested in certain special sets, typically rectan- 
gles of the form \\''^(a Jt bj), and for such sets the compu- 
tation is easy. 

The definition of M„ as the vector of componentwise 
maxima actually corresponds to regarding § as an ex- 
treme observation if %ij>v nJ (tj) for some j. More gener- 
ally, one may define £ to be an extreme value if 
£Ev r (A ) for some AC(0,1)'', in which case #(-4) has an 
interpretation as a measure of the clustering of such 
extremes. Note that the original definition of extremes 
corresponds to letting J 4-((0,r,)X...X(0,r i ,)y. Alter- 
nately one may consider taking yl-(?i,l)X...X(/ rf ,l) 
which corresponds to defining £ as an extreme observa- 
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tion if £ij>v,j(tj) for all j. Yet another choice is A- 
{t:2//>c>. 



4, Examples 

Wc conclude with two examples, both involving 
bivariate stationary sequences. 

Example 4.1 Let {17,,} be an iid sequence, and put 
£,i=7}„, and £,2=max{ij,-i, tj„ }. Let F denote the distri- 
bution of £,K£.i. £a) with marginals F { and F 2 . Then 
FiUWM&^.rWHifc ,Sjc, r,^x}~Ff(x) and 



v(cA)=c-"v(A) 



(15) 



F<jr t ,J 2 W {&,<*,, £.2 






if j:, ^jt 2 
if x,<Xi. 



If M'j) satisfies F,"(v„(r,))-«j, j'=l,2, then lim^« 
Fr(v„i(f 2 ))-f2 l,: so that v„ 2 (h)=Vn{il% Moreover 
f„i(/,)>v„2(f 2 ) if and only if t^t}' 2 , and so 



//„<v,</))=P{A/ n ==v„(0}->~ "(f)- 



•2 



The marginals of // are therefore //,(0=f| and f/ 2 (r 2 )= 
r 2 " 3 so that 0i-l and &- 1/2, and the dependence function 
of // is D H (t)-H{t u r 2 : )-f]Af 2 . For the associated iid 
sequence {£, > on the other hand, it is easily verified that 



Ut'2 , 11 'i<r 2 



1.2 

2 

1/2 



from which it follows that 






if/,>/ 2 " 2 

if r,<a 2 " 2 . 



Wc next consider a moving average sequence studied 
inRcf. [I6J. 

Example 4.2, Let {Zi-(ZnAj)'>, — «*:*«*, be a se- 
quence of iid random vectors in IR*. We assume the 
existence of a sequence of positive constants a,—^, and 
a measure v on IR 3 which is finite on sets of the form 
{x : |[x||>r }, r>0 (where }| • || denotes the Euclidean norm 
in IR 2 ), such that nP {a„ ' ZoG- H>* *>()- (Here W 
denotes vague convergence- of measures on IR 2 with 
respect to the metric rf(Xi,x 2 )=| ri" 1 — /v 1 Ivl 0i - 2 l, 
where for i- 1 ,2,r, and 0, denote the polar coordinates of 
X;, and a\/b=-max{a ,b }.) The measure v is necessarily 
of the form v({\ : H>r,0(x)EA})-r _ "S(A) for r>0 
and AC[0,27r), where 5(0 is a probability measure on 
10,2 tt) and a>0. Hence in particular [17], 



for all oO and all sets A with v(A)<x>. 

Define the bivariate moving average process 
\„~1J^CjZ„ ,, where {CjHcjjj]].-,}^ is a sequence of 
real 2X2 matrices satisfying SJ„ | c itt I *-c jc , £,/=! ,2, for 
some 5e<0,a). 8<1. For x=(.ri^ 2 )'GIR z , write A,/= 
{z: C,ze((-*vri)X(— =c r r 3 )) c }, where A' denotes the 
complement of a set A OR 2 . Then [16], 

lim P {a„ 'A/.<x}-exp{-y(x)}, and 

lim P {o;'A/„^xHexp{-7(x>}, xEIR 3 , 

where yW-V^viA^) and yW-KUjloA^). The ex- 
tremal index is therefore 0(x)=y(x)/f (x), xETR'. It fol- 
lows from the definition of A,j and Eq. (15) that this is 
in fact a function of XJX2. Note that the extremal index 
defined above differs from that in Sec. 3 in that it is 
defined on IR 2 rather than [0,1 1 J . However the two defi- 
nitions arc equivalent as may be seen by means of a 
suitable transformation from IR 1 to [0,1 J 2 . 

The actual calculation of 0(x) may be quite difficult 
in general, but possible to carry out under appropriate 
simplifying assumptions. 

Case (i). If Cy-c/T where C=[c u )tj_| and the tj-'s are 
non-negative constants, then A,j— cy'if(x) and 
U£A/-cfl(x), where B(x)={z : CzE((-*vr,)X 
(~ ac ^2>)' > and o.max{c, : fizQ}. Therefore by Eq. (15), 
v{A*j)-c?HBW) and KUjU^.j)-c"»'(B(*)) SO that 

Case (it). If the C/s arc diagonal, i.e., C/=diag[c,-i,Cja] 
with tj,>:0, /-1,2, then A,j-{z : Cjii^x, or c^z 2 >x 2 } 
and U^ SJ =((-=,T 1 /c 1 )X(-=o JC / r3 ))f where c,~ 
max{cy :f^0}, (-1,2. In particular, taking Xi-*^ and 
using Eq. (15) as in Case (i), we have p(A*j)- 
c",v((t : Zi>j:i}) ■ and v{V~ es A tJ }~c?v({z : Z\>*i}), so 
that the extremal index of {X ni } is ^-cfVi^c",. 

Case (Hi). Let D denote the support of v. If DC 
{2 : t,=0 or ^2=0} (which is the case if the coordinates 
of Zo are independent), then we may write 

L((-x^ 1 )x(-x t x 2 )f-a,xr n +fl 2 .r 3 ", Jc.^aO, (16) 

for suitable constants a^0 and a 2 ^0. Once again, as- 
suming the Cjjd's to be non-negative and writing 
cy-max{cjjj : j^0} for kj=l,2, we have (writing 
aAb=m'm{a,b}) 

A 1>; ni>-((-3c.r l /c ; , !1 A^ 2 /c J . 2 ,)X(-a:jr l /c,., 2 /a2/c;. 22 )) r n£) 
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and 

X(-3c^,/r„Aj:,/r ? ,>Vn'\ 
so that using Eq. (16) 



0(x)= 



ai(jf)/CuAJT; /c 2 i ) ""+a 2(jf l /c-i;Ay;/C2;) ^ 



Thus puttingjc 2 -^. wc have (^(o^n-Mit £>'(«! XJUt-ut 1 
oiSJoC/Si) and similarly, ^-(a,r^+^>' w)^"^"-^^^. ' 

If also c ; .[2-C/.2i-0 for each _/, (that i« if the f/r- ?te 
diagonal), then 



ew- 






and in particular, tfi—cfj/SJV/i t arid 8)*~t-wXJxfpi- Ml >^ 
that in this case the limiting distributions of M„ and A/„ 
are both independent, and hence (in accordance with 
Proposition 3.4) f?(x) is a convex combination of (,», sin.1 

The non-negative ncss of the C,*s assumed a'XM." is 
not crucial and may be relaxed, although it the c^" t f 
more involved calculations. 
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